Soft Hubbard gaps in disordered itinerant models with short-range interaction 
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We study the Anderson-Hubbard model in the Hartree-Fock approximation and the exact diag- 
onalization under the coexistence of short-range interaction and diagonal disorder. We show that 
there exist unconventional soft gaps, where the single-particle (SP) density of states (DOS) A fol- 
lows a scaling in energy E as A{E) oc exp[— (— 7log \E — Ef\Y] irrespective of electron filling and 
long-range order. Here, d is the spatial dimension, Ep the Fermi energy and 7 a non-universal 
constant. We propose a multi- valley energy landscape as their origin. Possible experiments to verify 
the present theory are proposed. 

PACS numbers: 71.27.+a, 71.23.-k 



Metal-insulator transitions (MIT) have been a funda- 
mental issue in condensed matter physics for a long time. 
The MIT is driven either by electron correlations, for ex- 
ample, as Mott transitions [H, or by random potentials 
as Anderson transitions Q . When the interaction causes 
an insulator, it opens a SP gap at the Fermi energy Ep- 
The Mott gap and a gap induced by an antiferromagnetic 
order (AF) are typical examples. On the other hand, 
the Anderson insulator exhibits absence of the gap with 
nonzero DOS at Ep-, where the insulators are character- 
ized not by the vanishing carrier number but by a singular 
relaxation time. This makes fundamental differences in 
low-energy excitations between the Anderson and Mott 
insulators. 

In real materials, however, electron correlations and 
randomness inevitably coexist, which may take on as- 
pects qualitatively different from the simple Anderson or 
Mott insulators 0, In particular, under the influ- 
ence of the interaction, Anderson insulators show quali- 
tatively different feature. Efros and Shklovskii [s^ (ES) 
have clarified that in the Anderson insulator with the 
long-range Coulomb interaction, a soft Coulomb gap 
opens in the SP DOS, A{E) with a power-law scaling 
as A{E) (X \E - Ef^ , a = d - I near Ep. The valid- 
ity of the ES theory was confirmed numerically and in 
experiments [sl later. In contrast, within the ES theory, 
short-range interactions do not generate soft gaps. 

Even for short-range interaction, however, soft gaps 
with a 0.5 were reported in a Hartree-Fock (HF) study 
in three dimensions (3D) [tI]. Recent numerical studies in 
two dimensions also show the suppression of DOS near 
Ep S @]- These suggest the presence of an unconven- 
tional mechanism which suppresses the DOS even with 
short-range interaction. In contrast, a numerical study 
with the dynamical mean-field theory (DMFT) claimed 
nonzero A{Ef) even in the insulating phases [lO|. Several 
mean- field studies gave similar results |lll,[l3l. We clearly 
need further studies for comprehensive understanding of 
the short-range case. 



Since the dielectric constant diverges at the MIT, ef- 
fects of the long-range part of the Coulomb interaction 
are restricted to low energies and the short-range part 
dominates electronic structures in the experimental en- 
ergy scale near the MIT. Therefore, unconventional soft 
gaps, if they exist, may be observed near MITs in real 
measurements. Indeed, recent photoemission results of 
SrRui-a^Tia^Oa [l3|, [ij] as well as a HF study with the 



long-range Coulomb interaction [15] indicate breakdown 
of the ES scaling near the MIT in 3D. 

In this letter, through numerical analyses of DOS at 
energies lower than those of the previous studies, we 
show even short-range interaction drives opening of a 
soft gap irrespective of the electron filling, originating 
from a mechanism entirely different from the ES the- 
ory. We call it soft Hubbard gap. We show numerical 
evidences of the soft Hubbard gaps with the HF ap- 
proximation in one and three dimensions, where DOS 
A{E) follows an unconventional scaling in energy E as 
A{E) oc exp[— (— 7log 1^ — Ef\Y] with 7 being a non- 
universal constant. Further support by the exact diag- 
onalization (ED) in one dimension (ID) is given. This 
scaling reduces to a power-law decay of A{E) toward Ef 
for d = 1 and even a faster decay for d > 1 in con- 
trast to the previous HF study [7]. To clarify the origin 
of the soft gap, we propose a phenomenological theory. 
The phenomenology is further numerically tested in de- 
tail against the HF results in ID. 

The Anderson-Hubbard Hamiltonian is defined by 



H 



on lattices with Ns sites and electrons, where t is a 
hopping integral, U the on-site repulsion, cj^ (Qct) the 
creation (annihilation) operator for an electron with spin 
a on the site i, riia = cj^c^cr and fi the chemical poten- 
tial. The random potential Vi is spatially uncorrelated 
and assumed to follow two models of the distribution 
Pv{Vi): the box type of width 2W, Pv{Vi) = 1/2W 
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FIG. 1: (color online), (a) Ground-state phase diagram of 
3D Anderson-Hubbard model at half filling for Gaussian dis- 
tribution of Pv. AFI, AF insulator; AFM, AF metal; PI, 
paramagnetic insulator (Anderson insulator); PM, param- 
agnetic metal, (b) DOS with system size 8 x 8 x 250: A 
(t = 1,U = 6,W = 5), B {t = 1,U = 4,W = 30), C 
(t = = 0,Ty = 30), D (t = 0,[/ = 4,1^ = 30). We em- 
ploy Lorentz broadening with a broadening factor 1.25 x 10""^ 
and 6.25 x 10~^ for A and B, respectively. The broken lines 
denote \E - Ef\ = 10"^ and 10"^ The DOS fits well with 
A(E) oc exp( — (— 7log \E — ^f|)^) shown by the fitting lines 
for 10~^ < 1^ — Ef\ < 10~^ as shown in the lower panel. 

{\Vi\ < W) with the average (Vi) = 0, and the Gaus- 
sian type, Pv{Vi) = ^exp{-V^^/2a^) {a^ = W^U). 
For both the distributions, fi = U/2 corresponds to half 
filling. We take the lattice spacing as the length unit. 

We first employ the HF approximation, where the wave 
function is approximated by a single Slater determinant 
consisting of a set of orthonormal SP orbit als {(/)n} is 
an orbital index). The HF equation reads: 



(2) 



where Ho is the one-body part of the Hamiltonian and 
we neglect (cJ|Q|). To find a site-dependent mean- field 
solution (riicr) for the HF equations, we employ the itera- 
tive scheme. One typically needs from several to several 
tens of initial guesses in obtaining convergent physical 
quantities such as AF order parameters and DOS. 

Figure [TJa) shows the ground-state phase diagram at 
half filling in 3D. We identify insulating phases by extrap- 
olation of the localization lengths ^ to the bulk limit. 
The localization length ^ is defined by the asymptotic 
behavior of SP orbitals near Ep at long distances as 
(j)n OC exp(— r/^), where r is the distance from the cen- 
ter of the orbital. We obtain the AF magnetic transition 
points by fitting the AF magnetic order parameter with 
the mean- field critical exponent, 1/2. Detailed analyses 
of the phase diagram will be discussed elsewhere. For 
U < 6 and with increasing W from 0, metals appear 
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FIG. 2: (color online), (a) DOS by HF in ID at t = 0.3, 
U = 1.0, W = 2.0, Ef = [7/2 - 1 (Ns = 14). Fitting of DOS 
gives a = 0.85 ± 0.07 (soUd line), which is in good agreement 
with the expected exponent of b' /b — 0.79 ± 0.02 (broken 
line), (b) Numerical estimates of P{R) and A{R). Fitting by 
Eqs. Jl and {|6| gives b' - 1.06 ± 0.01 and 6 = 1.34 ± 0.01. 
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FIG. 3: (color online), (a) DOS in ID with ED (open bound- 
ary condition): t = 0.1, U = 1.0, W = 1.0, Ef = U/2, 
(Ns = 2,3,4,5,6). We average the DOS over 3.2 x lO'^ 
realizations of disorder for Ns =6. (b) ScaUng plot by 
A{e,N-^) = N-^fieN^^"") {a = 0.075,0 = 0.375). 

from AFI as in the 2D result [16], with further reen- 
trant transition to insulators (AFI or PI). Naively one 
might expect A{Ey) > for W > 0. Figure [T](b) shows 
DOS for typical parameters. Indeed, there are no soft 
gaps when U or t is zero. However, we find soft Hub- 
bard gaps over the entire insulating phases in the case 
of U > and t > regardless of the AF magnetic or- 
der. Although a power law scaling A{E) ex \E — Ef\^ 
with exponents 0.5 < a < 1 looks fit in the range 
\E — Ef\ > 0.1 (not shown) being consistent with the 
previous HF study [7], closer look for \E — Ef\ < 0.1 fits 
better with A{E) oc exp(-(-7log \E-Ef\)^) with 7 > 
rather than the power-law scaling. 

The unconventional soft gaps exist also in ID regard- 
less of electron filling. Figure [2](a) shows DOS with the 
HF approximation for the box distribution of Py. Here, 
holes are partially doped with the chemical potential fi 
being shifted by —1.0 from the half filling to increase the 
average distance between electrons to capture long-range 
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(a) |0o) Ground state 

(b) |0i) Single-particle excited state 
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(c) |<^o) Nearly-degenerate state with |(/)o) 
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FIG. 4: (color online). Schematic illustration of (a) the 
ground state, (b) a SP excited state, (c) a nearly-degenerate 
state with the ground state and (d) a multiply-excited state, 
(e) Schematic of Vi dependence of excitation energies. 

asymptotic behavior easily. In contrast to the 3D case, 
they fit well with a power law A{E) (x\E — Ef\^ even at 
low energies. The gaps again vanish with the decreasing 
energy scale when t oi U becomes zero (not shown). 

In Fig. [Sja), we further show DOS with ED in 
ID. We assume a scaling function that A{e^N~'^) = 
N-^fieN^^"") = e^giN'^^'^e-^) corresponding to 
A{e,N-^ = 0) oc and A{e = O.N'^) ex N'^ (e = 
\E — Ef\)' As shown in Fig.[3](b), DOS well converges to 
this scaling function with a = 0.075 and P = 0.375. Al- 
though a possible logarithmic scaling cannot be excluded 
because of the small system size, the ED results are con- 
sistent with the HF results and support a mechanism of 
the soft gap working beyond the mean-field level. Be- 
cause the soft gap is restricted to very low energies in 
our ID study, further analyses at lower energies are de- 
sired in 2D, where only a pseudogap has been found so 
far 0. 

Now we discuss a possible origin of the soft gap. For 
simplicity without loss of generality, we restrict ourselves 
to a SP excitation for the electron side, namely, E > 
Ef' We consider the case of Hnt ^ where rint is the 
range of the interaction in the model. For t 7^ 0, virtual 
hopping of electrons generates effective interaction, which 
exponentially decreases with the mutual distance. This 
effect is not considered in the ES theory which regards 
electrons as classical particles. The DOS averaged over 
the random potential is obtained as 



A{E) 



Pv{Vi)Ai{E,Vi)dVi 



{Vt} 



(3) 



where the symbol {Vy} denotes a set of random potentials 
Vi except for Vi. Note that Ai{E, Vi) is the DOS under 
the condition of the fixed Vi at the site 1 and implicitly 
depends on {Vj}. Here we decompose the average over 
the random potential into the part for Vi as described 
by / iV(Vi)(iVi at fixed configurations of {Vj} and the 
subsequent average with respect to {Vj}. 

We discuss Fi-dependence of Ai{E^ Vi) for fixed {Vj}. 
When Vi decreases, the ground-state occupation of the 



site 1 changes from to 1 and then from 1 to 2 at Vid and 
Fic2, respectively. A possible ground state |(/)o) at Vi > 
Vici is illustrated in Fig. HJa), where the site 1 is empty 
and the total particle number is = Na and the energy 
Eo{Vi). Near Vici but for V > Vici, a SP excited state 
with Ne = Na -\- I and the energy Ei{Vi) is defined 
by the electron configuration except for the site 1 is fixed 
to be the same as |^o), as is illustrated in Fig.[Hb). One 
might think that becomes the ground state below 
Vici, where Ne = Na -\- I. In this case, however, the SP 
excitation gap Ei—Eq vanishes at Vici leading to absence 
of a gap in the Vi -averaged DOS at the site 1. Thus the 
numerical evidences of the soft gaps indicate that as 
a SP excited state is excluded by the electron correlation. 

In contrast to the ES theory, we assume a multi- valley 
energy landscape, which may be characteristic to random 
systems. Then there exist many arbitrarily- low-energy 
excited states whose configurations are the same with 
at the site 1 but globally different on other sites. In 
Fig. lU^c), we illustrate a state \(j)Q) at a local minimum 
Eq nearly degenerate with |(/)o), whose configurations not 
only at the occupied site n nearest to the site 1 at the dis- 
tance R but also farther sites (> R) are relaxed one after 
another. Figure [H^d) shows a SP excited state from 
1^0) with the energy E[ and the site-1 occupancy iden- 
tical with Here, the two nearly-degenerate states, 
and are separated by a barrier, where multi- 
particle relaxation is required to reach from one to the 
other. Now Ei is given by {Vi — Ep) + Uu + Eq^ 
where Uu is the interaction energy between electrons on 
the site 1 and those on the site i. Note that only the 
particles at the sites i which satisfy R < \i — 1\ ^ R^ i 
interact with the site 1 with the amplitude \Uii\ compa- 
rable to \Uin\ because of the localized nature. On the 
other hand, because £^0 — Eq and the configurations of 
on these sites are different from those of E[ is 
different from Ei by typically as much as |/7in|. Thus 
one can find 1^'^) with the energy E[ lower than Ei by 
as much as \Uin\ among many nearly-degenerate states 
with Now E[{Vi) and Eq{Vi) crosses at Vi = V-^^^ 
and for Vi < V-l^ the ground state becomes l^^)- Note 
that the excitation energy E[ — Ei is negative very near 
VYci but for Vi > V{^i. The state 1^'^) is not counted in 
DOS, because this state is not a SP excitation of |0o), 
but rather a multiply-excited state. Thus the energy dif- 
ference A = |£^i(^/ci) ~ ^i(^ici)l is lowest energy of 
SP excitations counted in A near Vi = V{^i. 

One might think that, as in the ES theory, it is possible 
to lower the energy of 1 0i ) from Ei to E[ by relaxing local 
electronic configurations only near the site n. It, how- 
ever, always increases the energy of the electrons other 
than those on the site 1, because they have already been 
optimized in the ground state and the increase dominates 
at large R. Thus a global reconstruction is required to 
lower the energy. 

From the above discussion, A scales as 
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A{R) = aexp{-bR), 



(4) 



where a and b are non-universal positive constants. Here- 
after we neglect logarithmic corrections. Under the as- 
sumption of linear dependence of the excitation energies 
on Vi as shown in Fig. [41(d), the local DOS averaged by 
Vi has a gap of A as follows; 

/CO 
Pv{Vi)A,{E,Vi)dVi (X H,{E-Ef-A), (5) 
viol 

where is the Heaviside step function. The same argu- 
ment applies around Vi = Vic2- 

The distribution of R with respect to {Vj} follows 

P{R) = a'exp(-6'i?^), (6) 

at long distances, where a\ h' are non-universal positive 
constants again. Equations (|4j) and (|6]) lead to 



Q(A) = P{Rm 



dR 



dA 



h' 



oc 



A-iexp(-^(-logA)'^), 



(7) 

where Q{A) are the distribution function of A. Equa- 
tions ^ and ((Tj) lead to 



A{E) oc 



/ 

J 



\E-Ef\ 



dAQ(A)«exp(-^(- 



\og\E-EF\f), 
(8) 

which is consistent with the observed scaling in ID and 
3D, further also in 2D (not shown) within the HE approx- 
imation. We also confirmed that this scaling is equally 
valid for a discrete distribution of Py (not shown). Eor 
d = 1, this leads to a power law with a no n- universal ex- 
ponent: A{E) oc \E — E-pf^^. Non-universal power-law 
distributions of energies without divergence of any length 
scales are common in Griffith phases [17]. Equation ^ 
indicates that a, namely the energy scale of the gaps 
vanishes as t or vanishes. Furthermore, the exponent 
a = b' /bis expected to decrease as t becomes smaller be- 
cause of the reduction of ^. These predictions are consis- 
tent with our HF results in ID. However, it conflicts with 
a DMFT study [13] and some mean-field studies 
which exhibit absence of the soft gaps. This 
because they ignore spatial correlations. The latter ig- 
nore inhomogeneity of the electronic structures. Indeed, 
a DMFT study with the intersite self-energy retrieves the 
suppression of DOS near Ep [9]. 

In Fig. [2](b), we show a further numerical evidence of 
our theory in ID. Figure [2](b) shows A{R) and P{R) cal- 
culated by the following procedure: First, we obtain the 
ground state for each realization of random potentials. 
We construct the lowest SP excited state by adding one 
electron to the lowest unoccupied orbital. Next we opti- 
mize the mean fields by the iterative scheme starting from 
those of the SP excited state with TVg fixed. Then A is ob- 
tained as the difference of these two excitation energies. 
We calculate R as the distance between the center of the 
lowest unoccupied orbital, r and those of the occupied 
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may be 



orbit als nearest to r in the ground state. We define the 
center of the orbital as the site which has the maximum 
weight. Fitting by Eqs. g]) and © gives b' = 1.06 ±0.01, 
b = 1.34 ±0.01. Estimated exponent ofV/b = 0.79 ±0.02 
is in good agreement with a = 0.85 ± 0.07 obtained di- 
rectly from DOS. This is a numerical evidence for the 
validity of our theory. 

Although the power-law was proposed to interpret the 
photoemission experiments 13|, [ij], our HF results in 
3D indicates that a different asymptotic behavior of the 
soft Hubbard gap emerges at lower energies, namely, 
< 10 meV. Since recent development of photoemission 
spectroscopy now allows us high-resolution measurement 
down to 1 meV, we believe that our paper provides incen- 
tive for such high-resolution photoemission experiments 
as well as for other measurement such as electrical trans- 
port measurement near the MITs. 

In summary, we have found an unconventional type 
of soft gaps in the Anderson-Hubbard model, although 
only short-range interaction is present. To clarify their 
possible origin, we have constructed a phenomenological 
theory. Detailed comparisons between our theory and 
the non-ES soft gaps observed in experiments are left as 
a future challenge. 
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